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(54) Method of measuring aberration in an optical imaging system 



(57) A method of determining aberration of an opti- 
cal imaging system comprises measuring at least one 
parameter, such as position of best-focus and/or lateral 
position, of an image formed by the imaging system. 



This is repeated for a plurality of different illumination 
settings of the imaging system, and from these meas- 
urements at least one coefficient, representative of ab- 
erration of said imaging system, is calculated. 
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Description 

[0001] The present invention relates to a method of measuring aberration in an optical imaging system, such as a 
Zgraphic projection apparatus. Observations are made as a function of parameters of the projection appara us, and 
5 thesTtKresence 'of different types of aberration can be quantified. The lithographic projection apparatus may 
comprise: 

a radiation svstem for supplying a projection beam of radiation; 
: a support structure for supporting patterning means, the patterning means serving to pattern the projection beam 
jo according to a desired pattern; 

. a substrate table for holding a substrate; and ....... 

. a projection system for projecting the patterned beam onto a target portion of the substrate. 

[0002] The patterning means here referred to should be broadly interpreted as referring to means that can be used 
to endow an incoming radiation beam with a patterned cross-section, corresponding to a pattern that is to .be created 
n a taVget portion of the substrate; the term "light valve" can also be used in this context Generally, the said pattern 
wii? correspond to a particular functional layer in a device being created in the target portion, such as an .ntegrated 
circuit or other device (see below). Examples of such patterning means include: 

- A mask The concept of a mask is well known in lithography, and it includes mask types such as binary, alternating 
phaL shift, and attenuated phase-shift, as well as various hybrid mask types. Placement of such a mask in the 
Nation beam causes selective transmission (in the case of a transmissive mask) or reflect™ in the case of a 
elective mask) of the radiation impinging on the mask, according to the pattern on the mask. In the case ofa 
mask the support structure will generally be a mask table, which ensures that the mask can be heW at a desired 
position in the incoming radiation beam, and that it can be moved relative to the beam if so desired. 
A programmable mirror array. An example of such a device is a matrix-addressable surface having a y.scoelast.c 
conSayeTand a reflective surface. The basic principle behind such an apparatus is that (for example) addressed 
areas of the reflective surface reflect incident light as diffracted light, whereas unaddressed areas reflect ^ncident 
St as undiffracted light. Using an appropriate filter, the said undirected light can be f ^ 
beam leaving only the diffracted light behind; in this manner, the beam becomes patterned according to the ad- 
Sing pa tern of the matrix-adressable surface. The required matrix addressing can be performed l us.n£ , su. ab e 
e Sronic means. More information on such mirror arrays can be gleaned, for example, from Un. ed States Paten s 
US 5 296 891 and US 5,523.193. which are incorporated herein by reference. In the case of a programmab.e 
mfrroSay the said support structure may be embodied as a frame or table, for example, which may be fixed or 

" . A programmable array. An example of such a construction is given in United States Patent US :5 229.87* 
which is incorporated herein by reference. As above, the support structure in th.s case may be embodied as a 
frame or table, for example, which may be fixed or movable as required. 
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For purposes of simplicity, the rest of this text may. at certain locations, specifically direct itself to examples involving 
a mask and mask table/however, the general principles discussed in such instances should be seen ,n the broader 
context of the patterning means as hereabove set forth. 

pool? Uhographic projection apparatus can be used, for example, in the manufacture of integrated c,cu,t .(ICs) 
" such a case the patterning means may generate a circuit pattern corresponding to an .nd.vidual layer of the IC and 
hi patirn Ln be imaged o'nto a target portion (e.g. comprising one or more dies) on a substrate s.hcon , wafeO that 
has been c 0 ^ t e d with a layer of radiation-sensitive material (resist). In general, a single wafer will contain a whole 
n!Lnrk of adjacent target portions that are successively irradiated via the projection system, one at a time. In curren 
ap^arls tfi£t££^ by a mask on a mask table, a distinction can be made between two different types of 
ES Sl^ Projection apparatus, each target portion is irradiated by exposing the en re mask 
rtte n o'n^he target portion in one go; such an apparatus is commonly referred to as a wafer stepper. In an a. ern Jve 
apparatus commonly referred to as a step-and-scan apparatus - each target port.on .s .rrad.ated by progressive y 
Sanntg the mask pattern under the projection beam in a given reference direction (the "scann.ng" direct.on) whrte 
onously scanning the substrate table paral.el or anti-parallel to this direction; since. ,n general, the projec .on 
system w5l have a magnification factor M (genera.ty < 1). the speed V at which the substrate table ,s scanned I w.N be 
a factor M times that at which the mask table is scanned. More information with regard to lithographic devices as here 
described can be gleaned, for example, from US 6.046.792. incorporated herein by reference. 
roOMI Ina manufacturing process using a lithographic projection apparatus, a pattern (e.g. in a mask) s imaged 
onto a substraTe Ta, is at least partia.iy covered by a layer of radiation-sensitive materia, (resist). Prior to th.s imag.ng 
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10 



15 



step the substrete may undergo various procedures, such as priming, resist coating and a soft bake. After exposure. 
Z substrate may be subjected to other procedures, such as a post-exposure bake (PES), development, a hard bake 
and measurement/inspeciion of the imaged features. This array of procedures is used as a bas.s to pafl em ar undMduri 
layer of a device e.g. an IC. Such a patterned layer may then undergo various processes such as etch ng. .on-.mplan- 
aL (doping), metallization, oxidation, chemo-mechanicat polishing, etc.. all intended to finish off an .ndmdual ^ayer 
rseveranayeVs are required, then the whole procedure, or a variant thereof, will have to be repeated for each new 
aver Eventually, an array of devices will be present on the substrate (wafer). These devices are then separated from 
one another by a technique such as dicing or sawing, whence the individual devices can be counted on a earner 
connected to pins. etc. Further information regarding such processes can be obta.ned. for example, from the book 
"Zochip Fabrication: A Practical Guide to Semiconductor Processing". Third Edition, by Peter van Zant. McGraw 

Hill Publishing Co.. 1997. ISBN 0-07-067250-4, incorporated herein by reference. 

moOS] For L sake of simplicity, the projection system may hereinafter be referred to as the "lens ; ha*»M» 
term should be broadly interpreted as encompassing various types of projection system, .ncluding refract.ve opt cs 
eflecfive optics, and catadioptric systems, for example. The radiation system may also include components operat ng 
according to any of these design types for directing, shaping or controlling the projection beam of rad.at.on, and such 
OTone 9 n2 ma'y a.so be referred to below, col.ective.y or singularly, as a "lens". Further, the , lithograph* apparatus 
maybeofatype having two or more substrate tables (and/or two or more mask tables). In such "mult.ple s^age dev,ces 
me additional tables may be used in parallel, or preparatory steps may be carried out on one or more tables wh.le one 
or more o°her tables are being used for exposures. Twin stage lithographic apparatus are descr.bed, for example. ,n 
20 US 5,969,441 and WO 98/40791. incorporated herein by reference. 

[0006] There is a desire to integrate an ever-increasing number of electronic components m an IC To realize this, ,t 
s necessary to increase the surface area of an IC and/or to decrease the size of the components. For the projection 
System this means that both the image field and/or the reso.ution must be increased, so that .ncreas.ng.y smaller 
STorlinewidths.can be imaged in a well-defined way inanincreasingly large image field. Th,s requires a projec .on 
25 system that must comply with very stringent quality requirements. Despite the great care with which such a Pr°Ject.on 
Stem s designed and me very high accuracy with which the system is manufactured, such a system may ^st. exh.b, 
Aberrations, such as sphericai aberration, coma and astigmatism. In practice, the projection system "lens") is thus not 
an Weal diffraction-limited system, but an aberration-limited system. The aberrations are dependent on pos,t.on n the 
Tmaoe field and are an important source of variations in the imaged line widths occurring across he .mage field, as 
30 Z I L influencing the focus, exposure latitude and so on. They a.so cause field-dependent overlay errors beNveen 
differen mask structures and/or different illumination settings. The influence of aberrations becomes .ncreas.ng.y s.g- 
£ wTth The application of newer techniques, such as phase-shift masks or off-axis illumination, to enhance the 
resolving power of a lithographic projection apparatus. 

rooOT] A further problem is that the aberrations are not constant in modern lithograph.c project.on systems^ In order 
35 to minimise low-order aberrations, such as distortion, curvature of field, astigmatism, coma and sphencal aberration, 
these pToIecSn systems generally comprise one or more movable elements. The wavelengths of the pro.ect.on beam 
or the pSL of the mask tab.e may be adjustable for the same purpose. When these adjusting fac.M.es are used 
o her smaller aberrations may be introduced. Moreover, since the intensity of the project.on beam must be s 
as possible the components of the projection system are subject to ageing so that the aberrafions may change durmg 
<o fhe Eme of the apparatus. Moreover/reversible changes, e.g. as caused by lens heating, may temporary change 

IwoT^Consequently there is a further problem of being able to measure the aberration reliably and accurately 
[0009] It is an object of the present invention to provide an improved method and apparatus for determ.n.ng aberrat.cn 

45 JuM^TcSwdinS.Te present invention provides a method of determining aberration of an optical imaging system 
comprising: 

- a radiation system for supplying a projection beam of radiation; . 

- a support structure for supporting patterning means, the patterning means serving to pattern the projection beam 
50 according to a desired pattern; 

a substrate table for holding a substrate; and 

- a projection system for projecting the patterned beam onto a target portion of the substrate, the method compns.ng 

the step of: 

55 patterning the projection beam with said patterning means; and 

characterized by the steps of: Hi« Q «»«t 
measuring at least one parameter of an image formed by the projection system, for a plurality of different 
settings of said radiation system and/or said projection system; and 
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calculating at least one coefficient, representative of aberration of said imaging system, on the basis of said 
at least one parameter measured at said plurality of settings. 

[001 1] Preferably, said plurality of different settings comprise different numerical aperture settings and/or sigma set- 
ings illumination modes or telecentricity modes; furthermore, one may use various types and sizes of test structures 
on for instance, one or more masks, to create different diffraction effects in the projection system. All such variation 
should be interpreted as falling within the meaning of the phrase "different illumination settings" as used in this text. 
The term "sigma (o) setting" refers to the radial extent of the intensity distribution in the beam at a pupil in the imaging 
svstem through which the radiation passes, normalized with respect to the maximum radius of the pupil. Thus, a sigma 
value of 1 represents an illumination intensity distribution with a radius at the pupil equal to the maximum radius of the 
nuDil The term "illumination mode" denotes the spatial distribution of the radiation at the pupil, which may be. for 
example disc-shaped, annular (which would be characterized by sigma inner and sigma outer settings), quadrupolar, 
dipolar soft-multipolar (including some radiation flux in between the poles), etc. The term "telecentricity modes en- 
compasses configuring the imaging system telecentrically and/or with varying degrees of non-telecentric.ty. for example 
by the use of prisms on top of a mask to tilt the illumination profile. These different settings can be selected conveniently 
in a lithographic projection apparatus. 

1001 21 The measured parameter can be one or more of: the position of best-focus of said image; the lateral position 
of said image; the deformation of said image; and other properties of exposing said image lithographically, such as 
line width and shape, and distance between adjacent structures. 
20 [001 3] Preferably, the plurality of different settings are selected such that the variation in the or each at least one 
measured parameter is substantially maximized. In this way the accuracy of the determined coefficient(s) can be im- 

raoiS' Preferably, the plurality of different settings are selected such that the variation in said at least one measured 
parameter resulting from aberration represented by one or more of said coefficients is substantially zero, whilst the 
variation in said at least one parameter as the function of a coefficient that is to be determined, is non-zero This 
technigue enables different aberration coefficients, such as Zernike coefficients, to be obtained independently of each 

[0015] The invention also provides a lithographic projection apparatus for projecting a patterned beam of radiation 
onto a substrate provided with a radiation-sensitive layer, the apparatus comprising: 

a radiation system for providing a projection beam of radiation; 

- a support structure for supporting patterning means, the patterning means serving to pattern the project.on beam 
according to a desired pattern; 
a substrate table for holding a substrate; 

- a projection system for projecting the patterned beam onto a target portion of the substrate; and 

illumination setting means for providing a plurality of different illumination settings of said radiation system 
and/or said projection system; 

characterized by further comprising: 
measuring means for measuring at least one parameter of a projected image formed by the projection system; 
control means for selecting a plurality of different illumination settings at which said measuring means takes 
measurements; and 

calculation means for calculating at least one coefficient, representative of aberration in said projection and/ 
or radiation system, on the basis of said at least one parameter measured by said measuring means. 

[0016] According to a further aspect of the invention there is provided a device manufacturing method comprising 
the steps of: 

(a) providing a substrate that is at least partially covered by a layer of radiation sensitive material; 
50 (b) providing a projection beam of radiation using a radiation system; 

(c) using patterning means to endow the projection beam with a pattern in its cross-section; 

(d) using a projection system to project the patterned beam of radiation onto a target portion of the layer of rad.at.on- 
sensitive material, 

55 and characterized by the steps of: 

measuring, prior to step (d), at least one parameter of an image formed by the projection system, for a plurality of 
different settings of said radiation system and/or said projection system; 
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calculating at least one coefficient, representative of aberration of said projection and/or radiation system, on the 
basts of said at least one parameter measured at said plurality of settings; 

correcting for said aberration on the basis of said at least one calculated coefficient, to reduce aberration of an 
image projected by said projection system. 

5 

[0017] Although specific reference may be made in this text to the use of the apparatus according to the invention 
in the manufacture of ICs, it should be explicitly understood that such an apparatus has many other possible applica- 
tions. For example, it may be employed in the manufacture of integrated optical systems, guidance and detection 
patterns for magnetic domain memories, liquid-crystal display panels, thin-film magnetic heads, etc. The skilled artisan 
10 will appreciate that, in the context of such alternative applications, any use of the terms "reticle", "wafer" or "die" in this 
text should be considered as being replaced by the more general terms "mask", "substrate" and "target portion", re- 
spectively. 

[001 8] In the present document, the terms "radiation" and "beam" are used to encompass all types of electromagnetic 
radiation, including ultraviolet radiation (e.g. with a wavelength of 365, 248, 193, 157 or 126 nm) and EUV (extreme 
is ultra-violet radiation, e.g. having a wavelength in the range 5-20 nm). 

[0019] Embodiments of the invention will now be described, by way of example only, with reference to the accom- 
panying schematic drawings in which: 

Figure 1 depicts a lithographic projection apparatus according to an embodiment of the invention; 
20 Figure 2 is a table summarizing the relationship between different low-order aberrations and their respective 

Zernike coefficients; 

Figures 3(a) and 3(b) are tables of computed best-focus positions (with respect to an aberration-free scenario) at 

different illumination settings as a function of aberration due to Zernike coefficients Z9 and Z16; 

Figure 4 is a plot showing the correlation between measured best-focus position and best-focus determined on 

25 the basis of Zernike aberration coefficients Z9 and Z16 measured according to the present invention, whereby 

each point in the plot corresponds to one of the illumination settings in Figure 3. The measured best-focus position 
is plotted along the vertical axis, in nm; the calculated best-focus position is plotted along the horizontal axis, in nm; 
Figures 5(a) and 5(b) are graphs showing the correlation between a Zernike coefficient measured using the method 
according to the present invention (plotted along the vertical axis, in nm) and as determined using another method 

30 (plotted along the horizontal axis, in nm), for Z9 and Z16 respectively; 

Figures 6(a) and 6(b) show plots of Zernike coefficients Z9 and Z16, respectively; as a function of x position along 
the slit of the projection system in a step-and-scan apparatus (in which the slit is scanned in the y-direction), as 
determined both by the method of the present invention and by another method. The vertical axis of both figures 
6(a) and 6(b) represents the value of the aberration coefficient, in nm. The horizontal axes represent the x position 

35 along the slit, in mm; 

Figure 7 is a graph of Z12 (shown along the vertical axis, in nm) as a function of x position along the slit of the 
projection system (shown along the horizontal axis, in mm) both as measured according to the method of the 
present invention and as measured using another method; 

Figures 8(a) and 8(b) are contour plots of best-focus position as a function of the outer sigma setting (shown along 
40 the horizontal axes) and numerical aperture (shown along the vertical axes) of the projection system as a result 

of aberration due to non-zero Zernike coefficients Z9 and Z16, respectively; 

Figure 9 is a contour plot of image contrast as a function of outer sigma setting (shown along the horizontal axis) 
and numeral aperture (shown along the vertical axis) of the projection system; 

Figure 10 is a table of calculated lateral shift in image position resulting from comatic (coma) aberration (non-zero 
45 Z7) as a function of illumination setting; 

Figure 11 is a graph showing correlation between measured x-shift (shown along the vertical axis, in nm) and 
calculated x-shift (shown along the horizontal axis, in nm) due to comatic aberration, whereby each point in the 
plot corresponds to one of the illumination settings in Figure 10; 

Figures 12(a) and 12(b) are contour plots of x-shift due to non-zero Zernike coefficients Z7 and Z14, respectively, 
so as a function of outer sigma setting (shown along the horizontal axes) and numerical aperture (shown along the 

vertical axes) of the imaging system; 

Figures 13(a) and 13(b) are plots of Zernike coma coefficients Z7 and Z8 (shown along the vertical axes, in nm), 
respectively, as functions of x position along the slit of the imaging system (shown along the horizontal axes, in 
mm), both as measured according to the method of the present invention and as measured using another method; 
55 and 

Figure 14 is a plot corresponding to that of Figure 13(a), but with a wavelength correction applied to the data of 
the present invention. 
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[0020] In the Figures, corresponding reference symbols indicate corresponding parts. 
First Embodiment 

5 [0021] Figure 1 schematically depicts a lithographic projection apparatus according to a particular embodiment of 
the invention. The apparatus comprises: 

a radiation system Ex, IL, for supplying a projection beam PB of radiation (e.g. UV or EUV radiation). In this par- 
ticular case, the radiation system also comprises a radiation source LA; 
10 a first object table (mask table) MT provided with a mask holder for holding a mask MA (e.g. a reticle), and connected 

to first positioning means for accurately positioning the mask with respect to item PL; 

a second object table (substrate table) WT provided with a substrate holder for holding a substrate W (e.g. a resist- 
coated silicon wafer), and connected to second positioning means for accurately positioning the substrate with 
respect to item PL; 

15 a projection system ("lens") PL (e.g. a refractive or catadioptric system or a mirror group) for imaging an irradiated 

portion of the mask MA onto a target portion C (e.g. comprising one or more dies) of the substrate W. 

As here depicted, the apparatus is of a transmissive type (i.e. has a transmissive mask). However, in general, it may 
also be of a reflective type, for example (with a reflective mask). Alternatively, the apparatus may employ another kind 

20 of patterning means, such as a programmable mirror array of a type as referred to above. 

[0022] The source LA (e.g. a Hg lamp or an excimer laser) produces a beam of radiation. This beam is fed into an 
illumination system (illuminator) IL, either directly or after having traversed conditioning means, such as a beam ex- 
pander Ex, for example. The illuminator IL may comprise adjusting means AM for setting the a-outer and o-inner values 
of the intensity distribution in the beam. In addition, it will generally comprise various other components, such as an 

25 integrator IN and a condenser CO. In this way, the beam PB impinging on the mask MA has a desired uniformity and 
intensity distribution in its cross-section. 

[0023] It should be noted with regard to Figure 1 that the source LA may be within the housing of the lithographic 
projection apparatus (as is often the case when the source LA is a mercury lamp, for example), but that it may also be 
remote from the lithographic projection apparatus, the radiation beam which it produces being led into the apparatus 
30 (e.g. with the aid of suitable directing mirrors); this latter scenario is often the case when the source LA is an excimer 
laser. The current invention and Claims encompass both of these scenarios. 

[0024] The beam PB subsequently intercepts the mask MA, which is held on a mask table MT. Having traversed the 
mask MA, the beam PB passes through the lens PL, which focuses the beam PB onto a target portion C of the substrate 
W. With the aid of the second positioning means (and interferometric measuring means IF), the substrate table WT 

35 can be moved accurately, e.g. so as to position different target portions C in the path of the beam PB. Similarly, the 
first positioning means can be used to accurately position the mask MA with respect to the path of the beam PB, e.g. 
after mechanical retrieval of the mask MA from a mask library, or during a scan. In general, movement of the object 
tables MT, WT will be realized with the aid of a long-stroke module (course positioning) and a short-stroke module (fine 
positioning), which are not explicitly depicted in Figure 1. However, in the case of a wafer stepper (as opposed to a 

40 step-and-scan apparatus) the mask table MT may just be connected to a short stroke actuator, or may be fixed. 
[0025] The depicted apparatus can be used in two different modes: 

1 . In step mode, the mask table MT is kept essentially stationary, and an entire mask image is projected in one go 
(i.e. a single "flash") onto a target portion C. The substrate table WT is then shifted in the x and/or y directions so 

45 that a different target portion C can be irradiated by the beam PB; 

2. In scan mode, essentially the same scenario applies, except that a given target portion C is not exposed in a 
single "flash". Instead, the mask table MT is movable in a given direction (the so-called "scan direction", e.g. the 
y direction) with a speed v, so that the projection beam PB is caused to scan over a mask image; concurrently, the 
substrate table WT is simultaneously moved in the same or opposite direction at a speed V= Mv, in which M is the 

50 magnification of the lens PL (typically, M = 1/4 or 1/5). In this manner, a relatively large target portion C can be 

exposed, without having to compromise on resolution. 

[0026] Aberrations of the projection lens in particular are considered in the embodiments of this invention. The pro- 
jection lens wavefront aberrations can be written as a series according to their angular form: 

55 

W(r t O)= £ R m (r)cos(m0)+R' m (r)sm(m0) 
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where rand 6 are radial and angular co-ordinates, respectively, (r is normalized) and m is an index indicating the 
contribution of the mth aberration. R and R'are functions of r. 
[0027] The aberration can also be expressed in terms of the Zernike expansion: 

W = Z/f,(r, 9) + Zjfj(r, G)+ Z/cf^(r,e)+... 

where each Z is the Zernike coefficient and each/is the corresponding Zernike polynomial. The function f take 
the form of the product of a polynomial in r and in sin or cos of me. For example, the comatic aberration can be 
represented by a Zernike series in Z7, Z8, Z14, Z1 5, Z23, Z24, Z34, Z35, etc., and, for example, the function associated 
with the Z7 coefficient [f 7 {r,Q ) in the notation above] is: 

(3r 3 - 2r)cos(6) 

[0028] The Zernike expansion for the lower-order aberrations is summarized in the table given in Figure 2. 
[0029] One effect of spherical aberration is to shift the position of best focus in the z direction and this is used in the 
present embodiment to determine Z9 and Z16. Firstly, the theoretical shift in best-focus position is calculated using 
known computational techniques, such as the program known as "Solid C", a commercial software package supplied 
by the company Sigma-C GmbH in Germany, for simulating and modelling optical lithography. Other suitable software 
packages such as one known as "Prolith", may alternatively be used. Figures 3(a) and 3(b) are tables giving the results 
of these calculations for two different lithographic projection apparatuses. The first two columns give the settings of 
the illumination mode in terms of the numerical aperture (NA), on the image side, and the sigma (a) setting. The 
illumination modes used in the tables of Figures 3(a) and 3(b) are annular illumination modes, and the two values in 
the sigma column give the inner and outer radii of each annulus; the difference between the two values gives the width 
of the annulus as a fraction of the maximum radius of the pupil. Sometimes the sigma inner value is quoted followed 
by the sigma outer value, and sometimes vice versa; the order of the sigma values is not significant. The range of 
radial positions of the beam at the pupil is related to the range of angular incidence of the projection beam onto the 
substrate 

[0030] The third and fourth columns in Figures 3(a) and 3(b) give the calculated best-focus positions in nm relative 
to a particular datum. The third column is calculated assuming that the lens has an aberration corresponding to a 
Zernike coefficient Z9 of 1 nm, the other Zernike coefficients being zero. The fourth column is for the same situation 
but with aberration due to Zernike coefficient Z16 of 1 nm, all other Zernike coefficients being zero. The values in the 
third and fourth columns (BF Z9=1nm and BF z16=1nm ) are effectively the gradients or partial differentials of the best-focus 
35 positions as a function of the respective Zernike coefficient, i.e.: 
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BF Z9 . iim (NA.cr)= — 



all other caeffidenti conj tan t 



in the approximation that 1nm is a small change in Z9, with a corresponding relation for Z16. 
[00311 The best-focus position is then measured using the actual lithographic projection apparatus at the same six 
different illumination settings as in each of the tables of Figure 3. The best-focus position is the z-position with maximum 
45 contrast e g defined by the maximum of a sixth-order polynomial fit to the contrast-versus-position curve as the position 
is moved from defocus, through focus and on to defocus. The best-focus can be determined experimentally using 
known techniques, such as the technique known as "FOCAL" (described below); alternatively, one may directly meas- 
ure the aerial image, e.g. using a Transmission Image Sensor (TIS) (described below), or an overlay microscope or 
commercial focus monitor. 

so [0032] FOCAL is an acronym for FOcus Calibration by using Alignment. It is a best-focus measurement technique 
for completely determining information about the focal plane using the alignment system of the lithographic apparatus. 
A special asymmetrically segmented alignment mark is imaged through focus on to a resist-coated wafer. The position 
of this imaged mark (latent or developed) can be measured by the said alignment system. Due to the asymmetric 
segmentation the position measured by the alignment system will depend on the defocus used during exposure, thus 

55 allowing determination of the best-focus position. By distributing these marks over the whole image field and using 
different orientation for the segmentation, the complete focal plane for several structure orientations can be measured. 
This technique is described in more detail in US 5,674,650, which is incorporated herein by reference. 
[0033] One or more transmission image sensor(s) (TIS) can be used to determine the lateral position and best focus 
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position (i e. vertical and horizontal position) of the projected image, under the projection lens. A transmission image 
sensor (TIS) is inset into a physical reference surface associated with the substrate table (WT). In a particular embod- 
iment two sensors are mounted on a fiducial plate mounted to the top surface of the substrate table (WT). at diagonally 
opposite positions outside the area covered by the wafer W. The fiducial plate is made of a highly stable material with 

5 a very low coefficient of thermal expansion, e.g. Invar, and has a flat reflective upper surface that may carry markers 
used with another fiducial in alignment processes. The TIS is used to determine directly the vertical (and horizontal) 
position of the aerial image, as projected by the projection lens, of a TIS pattern on the mask. It comprises apertures 
in the reflective surface close behind which is placed a photodetector sensitive to the radiation used for the exposure 
process To determine the position of the focal plane, the projection lens projects into space an image of a pattern 

10 provided on the mask MA and having contrasting light and dark regions. The substrate stage is then scanned horizon- 
tally (in one or preferably two directions) and vertically so that the aperture of the TIS passes through the space where 
the aerial image is expected to be. As the TIS aperture passes through the light and dark portions of the image of the 
TIS pattern the output of the photodetector will fluctuate (a Moire effect). The vertical level at which the rate of change 
of amplitude of the photodetector output is highest indicates the level at which the image of a TIS pattern has the 

15 greatest contrast and hence indicates the plane of optimum focus. The horizontal level at which the rate of change is 
highest indicates the aerial image's lateral position. An example of a TIS of this type is described in greater detail in 
US 4,540,277. Advantages of TIS include robustness and speed because it is a direct measurement technique not 
involving exposure of a resist. ( ... . . 

[0034] In cases where the present embodiment measures a parameter of the image (the best-focus position in this 

20 case) using the FOCAL technique at a plurality of different illumination settings, this will also be referred to by the 
acronym FAMIS, which stands for FOCAL At Multiple Illumination Settings. 

[0035] The actual Zernike coefficients of the projection lens system can then be extracted using the following relation, 
which gives the first-order correction to the best-focus position resulting from aberration due to the particular Zernike 
coefficients of interest: 



25 



BF meas (NA,c)= a + b-BF Z9= , nm {NA, o)+ c- SF Z16=1 nm (NA,o) 



30 



35 



40 



where BF meas is the empirical best-focus position measured as a function of NA and a. 
[0036] The coefficients b and c are called the Zernike coefficients Z9 and Z1 6, respectively; the coefficient a is related 
to Z4, and can be regarded as substantially constant in the current context. Consequently, the relation can be rewritten 
as the following equation: 

BF meas {NA,o) = const + Z9 BF Z9 ^ nm (NA, c)+Z16- BF 216s=1nm (/vAa) 

[0037] This equation is an approximation which neglects higher-order corrections, cross-terms involving both Z9 and 
Z16 and so on, although these could be taken into account if necessary. 

[0038] An equation, as above, is obtained for each illumination setting (NA and a values), six equations in total for 
the data in each of Figure 3(a) or 3(b). From these six simultaneous equations, three unknowns (const, Z9, Z16) are 
to be extracted As there are more equations than unknowns, the set of equations is over-determined. Multiple regres- 
sion or a least-square fit can be used to determine the best values of the Zernike coefficients. In matrix notation, the 
simultaneous equations can be written as: 



45 



50 



( BF ZMrm (NA l ,cr l ) BF^jNA,,*,) l) 
BF^jNA,,^) BF ZX6mX jNAt,<7 7 ) 1 



'Z9 ^ 

Z16 
\const j 



[0039] The first two columns of the matrix are the third and fourth columns of the tables in Figure 3(a) or 3(b) and 
55 are eigenvectors of the problem. 

[0040] In a more compact notation this can be rewritten as: 



8 



EP1 128 217 A2 



BF, 



= BF,-Z 



so Zernike coeffiecients are extracted according to: 



Z=BF " 1 BF n 



f0041] The Zernike coefficients Z9 and Z16 can be extracted in this way. Not only their (absolute) values are of 
nterest to the microlithographer. but also the effect they have on the best-focus position (and also on lateral shift or 
overlay error particularly in the case of other Zernike coefficients such as will be described with reference to the second 
embodiment) such that appropriate corrections can be made to the imaging system at different illumination settings. 
Sure 4 is a correlation plot between the actual measured best-focus position and the calculated best-focus position 
based on Z9 and Z16 determined according to the present embodiment of the invention. Each square represents a 
different illumination setting of the apparatus. As can be seen, there is excellent correlation, with the maximum focus 
the fitted line always less than 10 nm. The correlation coefficient R 2 = 0.972 (R 2 = 0 represents no 
correlation and R 2 = 1 represents perfect correlation). 

f0042] Further confirmation of the usefulness of the technique of the present invention is illustrated in Figures 5(a) 
and 5 b) These plot values of Z9 and Z16. respectively, as measured according to the present embodiment of the 
invent on (FAMIS) against values of Z9 independently measured by the lens manufacturer using a technique referred 
to hereafter as SIF (Shearing InterFerometry). SIF is a through-the-lens (TTL) interferometric technique developed by 
the lens manufacturer. Each square represents measurements at a specific x position along the slit of the imaging 
system in a step-and-scan apparatus, and each graph shows a best-fit straight line. The correlation is again excellent, 
with correlation coefficients R 2 of 0.86 and 0.83. The maximum deviation from the fit, which can be "lated to the 
p ecS of the measurement technique, is approximately 1 nm Zernike. This is well below ft. toe ranee wUh which 
the lens system is manufactured. The technique of the present invention has the advantage that it can be performed 
with the imaging system in situ in a lithographic projection apparatus. 

70043] Figure 6(a) and 6(b) show plots of Zernike coefficients Z9 and Z16 measured as a function of position ,n the 
x-direction i.e. perpendicular to the scan direction in an apparatus used in scan mode as described above. The aber- 
raJion oata we re measured on a step-and-scan system, but used in step mode. The aberration da a are pW to 
each x-position (along the slit), but averaged along the scan direction (y-direction). so the coefficients rep resent Mnte- 
grated values of Z9 and Z16. Each graph shows a plot of the Zernike coefficient measured according to the present 
FAMIS technique and as measured independently by the SIF technique. The graph in Figure 6(a) is a p\oi of the Z9 
coefficient and the graph in Figure 6(b) is a plot of the Z16 coefficient. These graphs again show the good correlat on 
oe JLTn the two measurement techniques. The correlation plots of Figures 5(a) and 5(b) correspond to the data ,n the 
oraDhs of Figure 6(a) and 6(b) respectively, with regard to said x-position. 

?0044] Figure 7 shows results of measurements of the Z12 Zernike coefficient using the FAMIS technique of he 
present invention and using SIF. Z12 (which relates to the presence of astigmatism: see Figure 2) can be calculated 
S "^present invention in the same way as Z9 or Z16, that is by computing the variation in astigmatism values 
(difference in best-focus position between x- and y- oriented lines) for a fixed amount of aberration due to, for example, 
a Z12 coefficient of 1 nm. The actual astigmatism values are measured at the same illumination settings, and a series 
of simultaneous equations formed and solved for Z12. . 
r00451 The FAMIS technique of the first embodiment of this invention can be used for extracting Zernike coefficients 
corresponding to even aberrations (i.e. aberrations with m = 0 or an even integer, see Figure 2). such as ^sphenca 
aberration as in Figures 3 to 6 and astigmatism as in Figure 7. because these aberrations influence the position of 
best-focus. Odd aberrations will be discussed hereafter in the Second Embodiment. 

[00461 Further refinements can be adopted to increase the accuracy of the measurement. Figures 8(a) and 8(b) are 
contour plots mapping the calculated best-focus position BF (from simulations) as a function of the NA and outer sigma 
va"ue oHhe illation setting, as a result of aberration represented by Zernike coefficients Z9 and Z16 
in Figures 8(a) and 8(b)). In each case the value of the respective Zernike coefficient was set at 10 nm. By choosing 
a range of different illumination settings that give rise to most change in best-focus position, the correspond.ng Zernike 
coefficient can be derived more accurately. . 
0047? For example, in Figure 8(a). if a series of measurements were taken with NA = 0.5 but at different outer ^sigma 
settings, very little variation in best-focus position would be measured, and consequently Z9 could only ^determined 
with relatively low accuracy. In contrast, a series of measurements at a fixed outer sigma setting of 0 7, bu t with £A 
Tanging from 0.5 to 0.7, would yield a larger range in best-focus position and more accurate determination of Z9. One 
technique would be to take a series of measurements at a range of illumination settings that ^J^^^ 
perpendicular to the contour lines in the plot of Figure 8(a). The same cons.derat.ons apply to determining Z16, and 
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clearly from the contour plot in Figure 8(b), a range of different illumination settings would optimize the determination 
of Z1 6. To optimize the determination of both 29 and Z1 6 simultaneously, the range of illumination settings should span 
variations in best-focus position as a function of both Z9 and Z16. 

[0048] The technique according to the present invention can also be enhanced if it is desired to determine particular 
5 Zernike coefficients independently of each other. By taking measurements at a series of illumination settings which 
follow, for example, a contour line in the plot of Figure 8(a), substantially no variation in best-focus position will be 
caused by Z9 aberration, and therefore substantially all the best-focus position variation will be due to the Z16 (or 
other) Zernike coefficients. In this way Z1 6 can be determined independently of Z9. The converse also applies, so that 
Z9 can be determined independently of Z16 by taking measurements at illumination settings corresponding to contour 
10 lines in Figure 8(b). 

[0049] A further consideration in the present embodiment is that for the FOCAL technique to measure the best-focus 
position with most accuracy, the contrast should preferably be as high as possible. Contrast is mapped in Figure 9 as 
a function of the NA and outer sigma values of the illumination settings. As can be seen in Figure 9, contrast is maximized 
at high NA settings and low outer sigma settings. 
is [0050] The parameter being measured can also be maximized by optimizing the structure being imaged, for example 
the type and size of the structure e.g. isolated lines and dense lines, isolated spaces, contact holes and so on, and by 
the choice of mask type, for example binary masks, attenuated masks and alternating phase-shift masks. 

Second Embodiment 

20 

[0051] The present invention can also be used to extract Zernike coefficients for odd aberrations. The present em- 
bodiment will be described in particular with reference to Zernike coefficient Z7, but could equally apply to other Zernike 
coefficients corresponding to aberrations with an odd value of m (see Figure 2). Odd aberrations result in a lateral shift 
of the image. 

25 [0052] Figure 1 0 is a table of computed lateral shift in the x-direction (also called x-shift, dX or AX) at eight different 
annular illumination settings; the settings of NA and sigma are given just as in Figure 3. The dX values are calculated 
on the basis of aberration due to a Zernike coefficient Z7 of 1 nm. 

[0053] Further calculations can give corresponding shifts due to the next-highest-order Zernike coefficient that results 
in x-shift, namely Z14 (see Figure 2). 
30 [0054] The lateral shift in image position can be measured in an actual lithographic projection apparatus using a 
distortion-measuring technique hereinafter referred to as "DISTO" (described below); alternatively, one may directly 
measure the position of an aerial image, e.g. using a Transmission Image Sensor (TIS) (described above) or other 
technique mentioned in connection with the first embodiment. 

[0055] DISTO is a technique for finding the lateral positional errors of the image plane (x,y plane). A standard align- 
35 ment mark, consisting of horizontal and vertical equal lines/spaces (e.g. line width 8 pm for the imaged mark), is imaged 
on to a resist-coated wafer. The position of this imaged mark (latent or developed) can be measured by the alignment 
system of the lithographic apparatus. The deviation with respect to the expected position is the distortion. By distributing 
these marks over the whole image field, the complete image field distortion can be measured. 
[0056] Where the second embodiment of this invention uses the DISTO technique at multiple illumination settings, 
40 this will be referred to as "DAMIS", an acronym for DISTO At Multiple Illumination Settings. 

[0057] The values of Z7 and Z14 can then be determined from the following relation, which gives the first-order 
correction to the measured shift in x-position resulting from X-coma due to non-zero Zernike coefficients Z7 and Z14: 

^meas(^^)= a+b ^ X Z7=lnm(W>A,a)+CdX 2l4=inm (/V>A,a) 

45 

[0058] The constants b and c in this case are now the Zernike coefficients Z7 and Z14, respectively, and can be 
obtained using the same matrix manipulation and multiple regression or least-square fitting as explained above in the 
first embodiment with reference to Z9 and Z16. In this case, term a is related to Z2, and is substantially constant. 
so [0059] Figure 11 shows a plot of correlation between experimentally measured X-shift (plotted along the vertical axis, 
in nm) and the calculated shift caused by a Zernike coefficient Z7 of 1 nm (plotted along the horizontal axis, in nm). 
The plotted points and calculated shifts correspond to the illumination settings in the table of Figure 10. The gradient 
of the line of best-fit gives the extracted value of Z7 for this apparatus, which in this particular example gives Z7 
approximately equal to 7 nm. 

55 [0060] Figures 12(a) and 12(b) show contour plots that map the shift dX in x-position in the illumination mode plane 
analogously to Figures 8(a) and 8(b). Here, the illumination setting was annular, with a ring width (sigma-outer) - (sigma- 
inner) = 0.3. Again, by appropriate choice of illumination settings at which to take measurement, the accuracy of de- 
termination of Z7 or Z14, respectively, can be improved, or variation due to one or other of Z7 and Z14 can be reduced 
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or eliminated to separate its contribution to the aberration from that of the other coefficient, as explained above for Z9 
and Z1 6 with respect to Figure 8. 

[0061] Figures 1 3(a) and 1 3(b) give further proof of the technique of the second embodiment of the invention. These 
Figures show plots ofZ7 (X-coma) and Z8 (Y-coma) as a function of x-position along the exposure slit. Values of the 
5 respective Zernike coefficients are plotted as obtained by the present invention and as independently measured using 
SIR The correlation is again good. 

[0062] The discrepancy between the values in Figure 13(a) can largely be reconciled by a tilt correction giving a 
variation in Z7 with x-position along the slit. Fig. 14 plots the same data as Figure 13(a), but with the DAM IS measure- 
ments of Z7 corrected by 0.39 nm per mm in x-position along the slit. This correction is achieved by an appropriately 
10 chosen wavelength shift. By way of illustration of the utility of the present invention, the total comatic aberration in the 
projection system indicated by the value of Z7 in Figure 1 3(a), measured according to the invention, could be corrected 
by a (slightly longer) wavelength shift. 

[0063] Although, in the embodiments described above, the sequence of operations consisted of calculating the ex- 
pected shift in best-focus position or in x- or y-position using a simulation for different illumination settings, followed by 

15 measuring the actual shifts at those illumination settings, it may of course be more practical to take the measurements 
first at particular illumination settings and then to calculate the appropriate quantities using the simulation at the actual 
illumination settings at which measurements were taken. Preferably, a sequence of test exposures is made sequentially, 
each with different NA and sigma settings, but with all other parameters of the exposures unchanged. 
[0064] To correct for the effects of machine drift, the measurement using the first illumination setting can be repeated 

20 at the end, and possibly also once or more often between the intervening different illumination settings, of a given 
measurement sequence. „.,-.. 
[0065] When making test exposures on wafers, it is advantageous to employ a technique known as micro-stepping 
or "die-in-die" exposures, because the same part of the wafer is used for all illumination settings, so all exposures are 
subject to substantially identical resist processing and metrology errors, and therefore the reproducibility of the test is 

25 increased. . . , 

[0066] The embodiments above have described lens aberration in terms of Zernike polynomials, but this is merely 
one of many possible sets of functions which can be used to describe wavefront aberrations, and the present invention 
can be used to measure aberrations in terms of parameters other than Zernike coefficients. 

[0067] Whilst specific embodiments of the invention have been described above, it will be appreciated that the in- 
30 vention may be practised otherwise than as described. The description is not intended to limit the invention. In particular 
it will be appreciated that the invention may be used for the projection system of a lithographic apparatus and in any 
other type of imaging system in which aberrations are to be measured independently of each other, reliably and accu- 
rately. 
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Claims 

1. A method of determining aberration of an optical imaging system comprising: 

40 - a radiation system for supplying a projection beam of radiation; 

- a support structure for supporting patterning means, the patterning means serving to pattern the projection 
beam according to a desired pattern; 

a substrate table for holding a substrate; and 

- a projection system for projecting the patterned beam onto a target portion of the substrate; the method com- 
45 prising the step of: 

patterning the projection beam with said patterning means, and 

characterized by the steps of: 
measuring at least one parameter of an image formed by the projection system, for a plurality of different 
so settings of said radiation system and/or said projection system; and 

calculating at least one coefficient, representative of aberration of said imaging system, on the basis of 
said at least one parameter measured at said plurality of settings. 

2. A method according to claim 1 , wherein said plurality of different settings comprise different numerical aperture 
55 settings. 

3. A method according to claim 1 or 2, wherein said plurality of different settings comprise different settings of pupil 
filling at a pupil plane in the imaging system. 
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4 A method according to claim 1, 2 or 3, wherein said plurality of different settings comprise different illumination 
modes selected from the group comprising disc-shaped, annular, quadrupolar, dipolar and soft-multipolar. 

5. A method according to any one of the preceding claims, wherein said plurality of different settings comprise tele- 
centric and non-telecentric illumination settings, or a plurality of non-telecentric illumination settings. 

6. A method according to any one of the preceding claims, wherein said at least one parameter is the position of best 
focus of said image. 

7. A method according to any one of the preceding claims, wherein said at least one parameter is the lateral position 
of said image. 

8. A method according to any one of the preceding claims, wherein said plurality of different settings are selected 
such that the variation in the or each at least one measured parameter is substantially maximized. 

9 A method according to any one of the preceding claims, wherein said plurality of different settings are selected 
such that the variation in said at least one measured parameter resulting from aberration represented by one or 
more of said coefficients is substantially zero, whilst the variation in said at least one parameter as the function of 
a desired coefficient is non-zero. 

10 A method according to any one of the preceding claims, further comprising the step of calculating a theoretical 
variation, for each of said plurality of different illumination settings, of the or each measured parameter as a function 
of a small change in said at least one coefficient. 

11. A method according to any one of the preceding claims, wherein said at least one coefficient is calculated by 
multiple regression or a least-square fit of a set of simultaneous equations. 

12. A method according to any one the preceding claims, wherein said at least one coefficient is a Zernike coefficient. 

13. A lithographic projection apparatus for projecting a patterned beam of radiation onto a substrate provided with a 
radiation-sensitive layer, the apparatus comprising: 

a radiation system for providing a projection beam of radiation; 

- a support structure for supporting patterning means, the patterning means serving to pattern the projection 
beam according to a desired pattern; 

a substrate table for holding a substrate; 

- a projection system for projecting the patterned beam onto a target portion of the substrate; and 

. illumination setting means for providing a plurality of different illumination settings of said radiation system 
and/or said projection system; 

characterized by further comprising: 

- measuring means for measuring at least one parameter of a projected image formed by the projection system; 

- control means for selecting a plurality of different illumination settings at which said measuring means takes 
measurements; and 

- calculation means for calculating at least one coefficient, representative of aberration of said projection and/ 
or radiation system, on the basis of said at least one parameter measured by said measuring means. 

14. An apparatus according to Claim 13, wherein the support structure comprises a mask table for holding a mask. 

15. An apparatus according to Claim 13, wherein the radiation system comprises a radiation source. 

16. A method of manufacturing a device using a lithographic projection apparatus, comprising the steps of: 

(a) providing a substrate that is at least partially covered by a layer of radiation-sensitive material; 

(b) providing a projection beam of radiation using a radiation system; 

(c) using patterning means to endow the projection beam with a pattern in its cross-section; 

(d) using a projection system to project the patterned beam of radiation onto a target portion of the layer of 
radiation-sensitive material, 
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and characterized by the steps of: 

measuring, prior to step (d). at least one parameter of an image formed by the projection system, for a plurality 
of different settings of said radiation system and/or said projection system; 

calculating at least one coefficient, representative of aberration of said projection and/or radiation system, on 
the basis of said at least one parameter measured at said plurality of settings; 

correcting for said aberration on the basis of said at least one calculated coefficient, to reduce aberrat.on of 
an image projected by said projection system. 

17 A method according to claim 16, wherein said correction step comprises correcting for changes in lateral and/or 
best-focus position of said image when changing between different illumination settings. 

18. A device manufactured according to the method of claim 16 or 17. 
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Fig.2. 
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Fig.5(a). *. 




Fig.5(b). 
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Fig.8(b) 
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Fig.10. 
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Fig.12(a) 
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Fig. 13(a). 20 
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